unclassified 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


la.  REPORT  SECURITY  CLASSIFICATION 

unclassified 


lb.  RESTRICTIVE  MARKINGS 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE 


3.  DISTRIBUTION/AVAILABILITY  OF  REPORT 


unlimited 


4.  PERFORMING  ORGANIZATION  REPORT  NUM8£R<S) 


MM  5488-86-19 


5.  MONITORING  f  GAN  I ZAT  ION  REPORT  NUMBER(S) 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

Mechanics  &  Materials  Ctr. 
Texas  A&M  University 


6b.  OFFICE  SYMBOL 
(If  applicable) 


7a.  NAME  OF  MC  ff  UNG  ORGANIZATION 

ONR  


6c.  AO  OR  ESS  (City,  State  and  ZIP  Code) 


College  Station,  Texas  77843 


7b.  ADDRESS  (City,  State  and  ZIP  Code) 


NAME  OF  FUNDING/SPONSORING 
ORGANIZATION 

ONR  


ftb.  OFFICE  SYMBOL 
(If  applicable) 


9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


Contract  N00014-86-K-Q298 


■c.  address  (City,  statt  and  zip  Co<u)  Mechanics  Division 
Office  of  Naval  Research/Code  1513A:DLU 
800  N.  Quincy  Street 
Arlington,  VA  22217  -  5000 


10.  SOURCE  OF  FUNDING  NOS. 


IT  TITLE  (Include  Security  Classification)  The .  D vnai" *  **  Energy 

Release  Rate  for  Two  Parallel  Steadily Propa¬ 
gating  Mode  III  Cracks  in  a  Viscoelastic  Body 


PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

element  no. 

NO. 

N° 

NO. 

4324  -520 

12.  PERSONAL  AUTHOR(S) 


13*.  TYPE  OF  REPORT 

Technical 

13b.  TIME  COVERED 

FROM  TO 

14  OATE  OF  REPORT  (Yr  .  Mo.,  Day) 

September  1986 

15.  PACE  COUNT 

29 

16.  SUPPLEMENTARY  NOTATION 

17. 


COS AT  I  COOES 


FIELD 

GROUP 

SUB.  GR. 

1 

18.  SUBJECT  TERMS  fConfmu*  on  reverie  if  neceuar y  and  identify  by  6loc*  number ) 

viscoelasticity  Barenblatt  zone 

dynamic  fracture  crack  growth 

dynamic  energy  release  rate 


19.  ABSTRACT  (Continue  on  reverie  if  neceieary  and  Identify  by  block  number ) 

The  dynamic  steady-state  propagation  of  two  parallel  semi-infinite,  mode  III 
cracks  is  considered  for  a  general,  infinite,  homogeneous,  and  isotropic  linearly  visco¬ 
elastic  body,  A  Barenblatt  type  -failure  zone  is  introduced  in  order  to  cancel  the  singular 
stress  and  a  formula  for  the  energy  release  rate  is  derived  which  provides  immediate 
comparisons  with  the  corresponding  single  crack  problem.  The  influence  of  crack  speed, 
crack  separation,  and  material  properties  upon  the  energy  release  rate  are  illustrated 
with  numerical  calculations  for  both  power-law  material  and  a  standard  linear  solid. 


20.  distribution/availability  of  abstract 

un class i f i e o/un limited  S  same  as  rpt.  □  otic  users  G 


21.  abstract  security  classification 

unclassified 


22*.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

Dr,  Richard  L.  Miller  ,  Code  1132P 


22b.  TELEPHONE  NUMBER 
(Include  Area  Code) 

(202)  696-4405 


22c.  OFFICE  SYMBOL 


DD  FORM  1473,  83  APR 


EDITION  OF  1  JAN  73  IS  OBSOLETE. 


unclassified 


SECURITY  CLASSIFICATION  OF  THIS  PAG 


The  Dynamic  Energy  Release  Rate  for  Two  Parallel 


Steadily  Propagating  Mode  III  Cracks  in  a  Viscoelastic  Body 


by 

L.  Schovanec*  and  J.  R.  Walton** 


Department  of  Mathematics 
Texas  A&M  University 
College  Station,  TX  77840  USA 


*  Permanent  Address:  Department  of  Mathematics,  Texas  Tech  University 
Lubbock,  TX  79409 

**  Supported  by  the  Office  of  Naval  Research  under  Contract  No. 
N00014-86-K-0298 


Mechanics  and  Materials  Center 
TEXAS  A&M  UNIVERSITY 
College  Station,  Texas 


library 


THE  DYNAMIC  ENERGY  RELEASE  RATE  FOR  TWO 
PARALLEL  STEADILY  PROPAGATING  MODE  III 
CRA  CKS  IN  A  VISCO  ELASTIC  BODY 


L,  SCHOVAN EC 
AND 

J.  R,  WALTON 


OFFICE  OF  NAVAL  RESEARCH 


DEPARTMENT  OF  THE  NAVY 


CONTRACT  N0001 4-86- K-0298 


WORK  UNIT  4324-520 


li  5  488-86  -19 


September  1986 


The  dynamic  steady-state  propagation  of  two  parallel  semi-infinite, 
mode  III  cracks  is  considered  for  a  general,  infinite,  homogeneous,  and 
isotropic  linearly  viscoelastic  body.  A  Barenblatt  type  failure  zone  is 
introduced  in  order  to  cancel  the  singular  stress  and  a  formula  for  the 
energy  release  rate  is  derived  which  provides  immediate  comparisons  with 
the  corresponding  single  crack  problem.  The  influence  of  crack  speed, 
crack  separation,  and  material  properties  upon  the  energy  release  rate  are 
illustrated  with  numerical  calculations  for  both  power— law  material  and  a 
standard  linear  solid. 


Introduction 


1  . 


In  [1]  the  problem  of  the  dynamic  steady-state  propagation  of  a 
semi-infinite  antiplane  shear  crack  in  a  linear  viscoelastic  body  was 
investigated.  By  considering  an  over-all  energy  balance  and  incorporating 
a  Barenblatt  type  failure  zone  at  the  crack  tip,  an  expression  for  the  rate 
of  energy  flux  into  the  crack  tip,  hereafter  referred  to  as  the  energy 
release  rate  (ERR),  was  derived  in  [2].  The  primary  objective  of  the 
latter  study  was  to  investigate  the  implications  of  using  the  ERR  as  a 
fracture  criterion  for  predicating  the  speed  of  dynamically  propagating 
cracks  in  a  viscoelastic  material.  In  particular,  a  steady-state,  dynamic 
elastic  calculation  would  predict  that  Mode  III  cracks  move  at  the  glassy 
shear  wave  speed.  However,  experimentally  observed  cracks  rarely  propagate 
faster  than  20%  of  that  speed.  Amplifying  on  the  elastic  calculation 
further,  what  one  finds  is  that  in  steady-state  the  stress  intensity  factor 
(SIF)  is  a  constant  independent  of  crack  speed,  while  the  ERR  is  a  monotone 
increasing  function  of  crack  speed  tending  to  infinity  at  the  shear  wave 
speed.  Thus,  increasing  crack  speed  corresponds  to  increasing  energy  flux 
to  the  crack  tip.  One  is  then  led  to  conclude  that  the  only  stable 
steady-state  speed  is  the  shear  wave  speed. 

For  viscoelastic  material  it  was  found  in  [1]  that  the  SIF  is  a 
monotone  decreasing  function  of  crack  speed  for  crack  speeds  between  the 
equilibrium  and  glassy  shear  wave  speeds  and  tends  to  zero  at  the  glassy 
wave  speed.  While  for  crack  speeds  less  than  the  equilibrium  wave  speed, 
the  SIF  has  the  constant  elastic  value.  In  [2]  it  was  shown  that  the  ERR 
has  much  more  complicated  behavior  and  in  general  does  not  depend 
monotonically  on  crack  speed.  In  particular  there  are  intervals  of  stable 
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and  unstable  crack  speeds  less  than  the  glassy  shear  wave  speed.  Thus 
viscous  effects  in  the  body  prevent  much  of  the  energy  input  by  the  applied 
loads  from  reaching  the  crack  tip  thereby  predicting  crack  speeds  much 
lower  than  results  from  an  elastic  calculation. 

It  has  been  conjectured  from  experimental  evidence*  that  another 
mechanism  limiting  crack  speeds  in  viscoelastic  material  is  the  interaction 
of  microcracks  around  the  tip  of  a  macrocrack.  As  a  first  step  toward 
investigating  this  conjecture,  the  problem  of  two  parallel,  semi-infinite 
Mode  III  cracks  dynamically  propagating  in  steady-state  is  considered  in 
this  paper.  As  demonstrated  below  it  is  found  that  the  energy  input  into 
the  tip  of  one  crack  may  be  reduced  by  up  to  a  factor  of  one-half  due  to 
the  prescence  of  another  crack.  The  influence  of  crack  interaction  is 
greatly  influenced  by  crack  speed,  crack  separation,  and  material 
properties  (e.g.  stress  relaxation  time  and  the  ratio  of  glassy  to 
equilibrium  shear  wave  speeds). 

In  the  next  section  the  relevant  boundary  value  problem  is  presented 
and  solved  by  reduction  to  a  Riemann-Hilbert  boundary  value  problem  in  a 
manner  similar  to  [1].  Section  3  contains  the  calculation  of  the  ERR 
following  the  method  of  [2].  The  paper  concludes  with  numerical  examples 
in  Sections  4  that  illustrate  the  general  results  of  Section  2  and  3. 


*  Private  communication  with  Professor  Wolfgang  Knauss. 
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2 .  The  Boundary  Value  Problem 


The  specific  boundary  value  problem  to  be  investigated  is  that  of  the 
steady  propagation  (to  the  right)  with  speed  V  of  two  parallel, 

^  semi-infinite,  anti-plane  shear  cracks  in  a  general,  homogeneous,  and 

isotropic,  linearly  viscoelastic  body.  The  shear  modulus,  y(t),  is 
assumed  only  to  be  a  positive,  non— increasing,  and  convex  function  of 
•  time,  t.  The  governing  equation  of  motion  for  the  out  of  plane 

displacement,  u^ ,  is 


*  dAu3  =  Pu3,tt 


where 
y  *  de 


A  is  the  two-dimensional  Laplacian,  A  =  (3^/3x^)  +  (S^/Sx^),  and 
denotes  the  Riemann-Stieltzes  convolution 


V  *  de  "  IL  M(t  -  T)de(t). 


The  cracks  are  assumed  to  lie  in  the  horizontal  planes  corresponding  to 
=  h  and  =  -h  with  their  tips  aligned  in  the  plane  x^  =  0.  If 

the  crack  faces  are  subjected  to  equal  loadings,  the  problem  is 
anti-symmetric  about  the  plane  y^  =  0.  With  the  introduction  of  the 
Galilean  variables,  x  =  x^  -  Vt,  y  =  x^,  it  then  suffices  to  solve  the 
equation  of  motion  in  the  half-plane  y  >  0  subject  to  the  boundary 
conditions 


°23(x’h+)  =  a23(x’h_)>  -“  <  x  <  00  (1) 

a23(x,h)  =  gl  (y  *  du3)  =  f (x ) ,  x  <  0 


(2) 


(3) 


u^(x,h4')  =  u^(x,h  ),  x  >  0 


u~(x,0)  —  0 ,  <  x  < 


(4) 


The  boundary  value  problem  is  solved  by  proceeding  as  in  [1].  The 

Fourier  transform,  defined  bv 


*  00  ipx 

f(p»y)  =  /  oo  e  f(x,y)dx 


is  applied  to  the  equation  of  motion  resulting  in  the  ordinary  differential 
equation 


u3(p,y) 


=  ^  2(p)u3 


2  2 

where  Y  (p)  =  p  +  ipVp/p(-Vp).  The  solution  of  the  differential  equation 
is  written  as 


A(p)e  y  >  h 

B(p )sinh( Y  y)  0  <  y  <  h 

2 

so  that  (4)  is  satisfied  and  a  square  root  of  Y  with  positive  real  part 
is  selected  to  insure  the  stresses  and  displacements  vanish  as  y  +  00 . 

From  (1)  it  easily  follows  that 

A(p)  =  -B(p)e^cosh(Yh) . 

Define  u(x)  =  u3(x,h+)  -  u3(x,h~)  and  a(x)  =  a23(x,h”)  =  a23(x,h+). 

Then  a  simple  calculation  gives 


u3(p>y)  =  y 
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°(p)  =  -( iVp)u(-Vp)Y(p)e  Yhcosh(Yh)u( p) . 


(5) 


If  the  Carson  transform  g(s)  of  the  function  g(t)  is  defined  by 


g(s)  =  g( 0)  +  /”  e  tSdg(t) 


and  Yj(p)  -  [1  -  PV2/p ( iVp) ] ^2,  (5)  may  be  written  as 


a(p)  =  -isgn(p)p(iVp)Y1(p)  (1  +  e_2h|  PIT  1  )u,  ^p) .  (6) 


Define 


G^p)  =  -isgn(p)u(iVp)Y1(p)  =  sgn(p)G  (p). 


G,(P)  4(1  +  e'2h  P  Yl)  -l6 


2  2 


(p) , 


and  G  =  G^.  From  (2)  and  (3)  it  then  follows  that  (6)  may  be  viewed  as 
the  Riemann-Hilber t  problem. 


•  F+(p)  =  G(p)F  (p)  +g(p)  (7) 

where  F+(p)  =  o+,  F“(p)  =  u“  g(p)  =  -a  =  -f(p),  and  f+(x)(f(x)) 

denotes  the  restriction  of  f(x)  to  x  >  0  (x  <  0).  Note  that  as  h  0, 

G -»•  G^  and  (7)  is  precisely  the  problem  of  a  single  crack  solved  in  [1]. 

Since  the  solution  of  (7)  is  constructed  in  a  similar  manner  to  [1],  only  a 

brief  outline  of  the  method  of  solution  is  presented  here. 

To  solve  (7)  first  consider  the  homogeneous  problem  of  finding 
±  \ 

functions  X  (z)  analytic  for  lm(z)  '  0  which  satisfy  the  boundary 


relation 
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X+(p)  =  (p)X-(p). 


(8) 


-f  -f  -f 

The  solution  of  (8)  may  be  written  as  X-  =  X~X~  where  x+(p)  = 
i(p)X"(p),  i  =  1,2.  X~  (  z  )  is  given  in  [1]  and  X*  =  w±X±  where 


X  (z)  =  exp  [r±(z>] 


+  i  ,co  i°e(g  9(t  ) ) 

r~(z)  =  ~  f  - - - 

2iri  1  t  -z 


dt  , 


(9) 


and  w  (p)/w  (p)  -  1/2.  Though  X_(z)  is  not  uniquely  determined,  (7) 


has 


a  unique  solution  given  by 


F  (z)  = 


X  (z)X  (z)  g(T )/(xtxt)(x  ) 


12' 


2tt  i 


t  -z 


di 


(10) 


Utilizing  (10)  and  some  trivial  modifications  of  the  methods  of  [1]  it  can 
be  shown  that  the  asymptotic  expansion  of  for  x  near  zero  is 

given  by  a +3(x)  -  K//x  where  the  stress  intensity  factor  K  is  defined  by 


K  = 


jG1  C00) 


1/2 


/if 


if" 

2tt  j 


g(x  ) 


x|(t)X+(t  ) 


dx 


(ID 


In  the  next  section  attention  is  focused  upon  the  calculation  of  the 
energy  release  rate,  the  primary  interest  of  this  paper.  For  special  forms 
of  the  crack  loadings  explicit  formulas  for  the  ERR  and  the  SIF  are 
developed.  The  paper  concludes  by  presenting  numerical  results  for  the 
special  cases  of  a  power-law  material  and  a  standard  linear  solid. 
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3 .  The  Energy  Release  Rate 

As  discussed  in  [2] ,  the  energy  flux  into  the  tip  of  an  extending 
crack  in  a  viscoelastic  material  depends  upon  the  entire  history  of  the 
singular  asymptotic  stress  and  strain  fields  at  the  crack  tip  during  the 
time  the  tip  is  advancing.  This  energy  flux,  referred  to  previously  as  the 
energy  release  rate  and  hereafter  denoted  by  G,  is  difficult  to  calculate 
from  the  singular  fields.  For  this  reason  it  is  convenient  to  introduce  a 
Barenblatt  type  failure  zone  behind  the  crack  tip  in  order  to  cancel  the 
singular  asymptotic  fields  in  front  of  the  crack.  Adopting  the  same 
notation  as  in  [2]  and  [3],  the  appied  external  tractions,  o~^(x,0)  =  f(x), 
will  now  be  denoted  by  o— (x)  with  the  cohesive  failure  stresses  o— (x) 
acting  in  a  failure  zone  of  length  af.  As  pointed  out  in  [2]  what  results 
from  incorporating  a  failure  zone  behind  the  crack  tip  is  that  the  energy 
flux  into  the  crack  tip  may  be  expressed  as 


G  =  -  /°  a~(x)u~,  (x)  dx 

3^  r  x 


(12) 


where  now  u(x)  =  u^Cx^"1")  -  u^(x,h“)  is  the  relative  crack  face 

displacement  corresponding  to  the  combined  loading  a~  +  a~ 

As  in  the  previous  studies,  the  external  load  a-  and  failure  zone 

e 

stresses,  will  be  assumed  to  have  the  forms 


x/ a  x/a, 

=  V  »  afM  "  "Lfe  ' 


<  x  <  0. 


(13) 


For  af/ae  <<  1  the  essential  features  of  the  Barenblatt  model  are 
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satisfied  by  the  assumptions  (13),  namely  a  set  of  cohesive  stresses  and 
associated  length  scale  a^  and  a  length  scale  a^  associated  with  the 
applied  load  such  that  cancels  the  singular  stresses  produced  by 

c~.  When  0“  and  o~  are  given  by  (13),  (12)  is  replaced  by 


=  -/: 


°f (x)u' 


(x)dx. 


(14) 


Applying  the  Parseval  formula  for  the  Fourier  transform,  (14)  may  be 
written  as 


G  =  -  f_m  o~(p)u  >x(p)dp 


05) 


where  f  denotes  the  inverse  Fourier  transform 


f(p)  =  27  /_ *  f(x)e  1XPdx. 


From  (13)  it  follows  that 


-afLf 

°f(p)  "  (1+ia.p)  ’ 


afLf 


°f(p)  2Tr(l-ia,p) 


06) 


with  similar  expressions  for  o  and  a  .  Since  u~  has  an  analytic 

e  e  9  x 

extension  to  the  lower  half-plane  given  by  F~(z)  and  has  a 

meroraorphic  extension  to  the  lower  half-plane  with  a  simple  polh  at  -i/a^, 
the  integral  in  (15)  may  be  evaluated  using  residues.  This  calculation 
yields 


G  =  LfF-(-i/af). 

Utilizing  the  general  method  for  evaluating  F"(-i/a^)  presented  in  [3] 


one  finds 
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G  = 


(af  -  ag)  L2  X~(-i/af )X+(i/af ) 


(a  +  a  ) 
t  e 


[X+(i/a  )] 
e 


(17) 


where  X  are  defined  by  (8).  A  convenient  factorization  of  G  which 

allows  for  comparisons  between  the  single  and  two  crack  problems  is 

+  +  + 

achieved  by  recalling  X  =  X^2  and  writing  G  =  G ^  where 


Uf  "  3e)  Le  ^(-i/af)X|(i/af) 
G1  (a,  +  a  )  2 


[X|(i/ae)]‘ 


and 


G2  " 


X;(- i/af)X+(i/af )  2X0(-i/af )X+(i/af ) 


f y  2' 


[X+(i/ae)]' 


[X+(i/ag)]‘ 


(18) 


It  was  shown  in  [1]  that  two  cases  arise  naturally  in  constructing 

X*:  0  <  V  <  C*  and  C*  <  V  <  C  where  C*  =  (p(“)/p)1/2  and 

1/2 

C  =  (U(0)/P)  are  the  elastic  shear  wave  speeds  corresponding  to  the 

equilibrium  and  glassy  values  of  the  shear  modulus  y(t).  If  a 
characteristic  time  scale,  is  introduced  and  the  shear  modulus  given 

in  the  form 

P(t)  =  M(t/x) 


where  ra(s)  is  a  nondimensional  function  of  s  with  lim  m(s)  =  1  and 

g->O0 

then  Y,(p)  =  (1  -  PV^/p^mC  iVT  p) )  *  ^  and  from  [2]  it  then 

£-KX)  i 


follows  that  may  be  expressed  as 
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IF  3  L 

=  (ilLj  (_£_£.)[ 

'  1+e J  K  "Jn  J  L 


2u„ 


1  ~  (£8) 


—  av  J 


2  1/2 


i-y2«21> 


(19) 


In  (19)  Che  following  nondimensional  parameters  have  been  introduced 


“  =  C*T/ae,  8  =  q0ae,  Y  =  V/C*,  £  =  af/a 


f  e 


and  qQ  =  0  for  0  _<  V  <  C  while  for  C  <  V  <  C ,  qn  satisfies 


a,  -qnvt  y 

Vq 0  J Q  (U(t)/u(0))e  dt  =  (V/C)  . 


Now  write  G 2  =  2G2  so  that  from  (18)  and  (9)  there  results 

los(V  -  ik  ‘  -  ^Jds. 


Making  the  substitution  s  =  p/ap  and  introducing  the  nondimensional 


parameter  h/a  =  one  finds 
e 


log  (G2) 


-  4  /liog[i+cOT(-2lpl«(i-r2/;(if£))1/2][J-  +  _/JV2-»  ]dp 


p2+l  (p2+l)(l+(£p)2) 


Vo  h<P>d» 


(20) 


where 


h(p) 


tan  ^ (v/u)p(£ ^-1 )  _  log  /u2  + 
( 1+P~)( l+(ep)^)  p^+l 


and 


1  1 


,  .  x  x 

U  =  1  +  e  cos  y,  v  =  e  sin  y 


X  =  -/2  p6(/s2+t2  +  s)1/2,  y  =  _/2  p6 (/ s 2+t  2 

(1  -  y2/^(HLP))1/2  =  (S  +  it)1/2. 


-  s) 


1/2 


(21) 


Note  that  the  effects  of  the  crack  separation  distance,  h,  are 
contained  in  the  term  G^  while  ,  which  was  computed  in  [2]  , 
represents  one  half  the  energy  flux  into  the  crack  tip  for  the  single  crack 
problem.  It  is  easy  to  verify  that  as  h  -*■  »,  G  2G  ,  the  ERR  for  the 
single  crack  problem,  and  as  h  +  0,  the  work  input  into  the  crack  tip  is 
reduced  by  one  half  since  G  G^  . 

The  methods  used  to  evaluate  (14)  when  applied  to  (11)  show  that  the 
stress  intensity  factor  may  be  expressed  as 


K  = 


-Le|V~> 


1/2 


/  7T 


x!(i/ae)  X+(i/ae) 


=  K 


1  Xt(i/a  ) 

z  e 


where  is  the  SIF  for  the  single  crack  problem  found  in  [1].  The  same 


reasoning  leading  to  (20)  yields 


log  (X+(i/ag))  =  j  fQ  k(p)dp  (21) 

where 


1  7  9  —  1 

k(p)  =  - (log  / u  +v  +  p  tan  1(v/u)) 

(1+P2) 

with  u  and  v  defined  as  above.  Again  it  is  easily  verified  that  as 
h  -*■  ®,  K  -*■  Kl  . 


12 


This  section  concludes  with  the  observation  that  the  corresponding 
results  for  the  elastic  problem  are  easily  recovered  from  the  above 
calculations  by  setting  a  =  0  and  n  =  0,  in  which  case  C  =  C*  and 
0  <  Y  <  1.  In  the  next  section  the  qualitative  behavior  of  the  ERR  is 
illustrated  by  numerically  evaluating  G  for  special  cases  of  the  shear 
modulus . 
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4 .  Numerical  Examples 


G  will  now  be  computed  for  the  special  cases  of  a  standard  linear 
solid  and  power-law  material.  First  considered  is  the  standard  linear 
solid,  which  is  modeled  by  a  shear  modulus  of  the  form 

U(t)  =  u^C  1  +  n  e"t/T  ) 

=  p^mCt/t  ). 

In  this  case  m(w)  =  ( l+w(  l-4n  ) )/( 1+w)  and  s  and  t  defined  in  (21)  are 
given  by 


s  = 


ilzl  )  +  (p°or)z(  l^n )( Hn-y2) 

1  +  (pay ( 1+p) ) 2 


t  = 


3 

pay  n 


1  +  (pay(l-m  ))' 


For  the  standard  linear  solid  restriction  0  <  V  <  C  corresponds  to 
0  <  y2  <  1-ki  . 

In  the  case  of  a  power  law  material  the  shear  modulus  is  assumed  to 
have  the  form 


u(t)  =  ujl  +  (t/r  )  n),  0  <  n  <  1 


=  y  rn(t/T  )  . 


Now  m(w)  =  1  +  r(l-n)wn  and  if  =  a [T ( 1-n) ]  1/n,  s  and  t  are  given 


by 
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1  +  (pct.T)n  cos(— “) 

s  =  1  -  y2{ - l- - 2 - -} 

1  +  ZCpo^y)11  cos(—“)  +  (pa  ^y)  2n 

Y  (p«jY)  sin  (— ) 

1  +  2(pot  )  cos  (  — )  +  (pot^Y) 

In  this  case  the  glassy  wave  speed  C  is  infinite  so  that  0  <  Y  <  00 . 

Since  a  formula  for  is  given  in  [2],  to  calculate  G  it  suffices 

to  numerically  integrate  (20).  Certain  care  must  be  exercised  in 
performing  the  necessary  integrations  as  the  integrand  h(p)  may  display  a 
steep  cusp  near  the  origin  along  with  severe  oscillation  for  certain  ranges 
of  the  parameters  considered.  Adaptive  Gaussian  integration  techniques 
were  utilized  in  obtaining  the  following  results. 

In  Figure  1,  G^  for  the  standard  linear  solid  is  plotted  against 
log  Y  when  Y  <  1  and  against  log((/q+l  -  l)/(/n+l  -  Y))  for 
1  <  Y  <  / 1+n  with  6=1,  n  =  10,  e  =  .01,  and  a  =  .1,  1,  10,  100. 

One  explanation  for  these  results  is  suggested  by  noting  that  6/a  =  h/c*r 
so  that  small  values  of  a  correspond  to  a  crack  separation  distance  that 
is  large  compared  to  the  distance  a  signal  travels  at  the  speed  of  an 
equilibirum  shear  wave  during  a  characteristic  unit  of  time.  Figure  1 
illustrates  that  when  6 /a  is  large,  the  power  input  into  the  crack  tip 
approaches  that  corresponding  to  the  single  crack  problem.  As  6 /a 
becomes  small,  this  effect  occurs  only  at  crack  propagation  speeds  that 
approach  the  glassy  shear  wave  speed.  The  normalized  ERR  for  the  standard 
linear  solid,  g,  defined  by 


1-e 

1+e 


A- 


a  L2 
e  e 


2y 


-Jg 


is  graphed  in  Figure  2.  The  behavior  of  g  is  qualitatively  similar  to 
that  for  the  single  crack  problem  and  the  appropriate  comments  in  [2]  apply 
here.  The  corresponding  results  for  the  power  law  model  appear  in  Figures 
3  and  4.  The  only  significant  difference  is  that  in  Figure  3  G always 
attains  the  value  corresponding  to  the  single  crack  problem.  This  is 
explained  by  the  fact  that  for  the  power-law  model,  the  glassy  shear  wave 
speed  is  infinite.  Thus  this  model  allows  sufficiently  large  crack 
propagation  speeds  to  overcome  the  effect  of  d/a,  which  may  be  regarded  as 
a  measure  of  the  influence  that  one  crack  tip  has  on  the  other. 

The  remaining  figures  apply  to  a  standard  linear  solid.  (The 
corresponding  graphs  for  the  power-law  model  are  qualitatively  similar.) 

In  Figures  5  and  6  is  graphed  versus  log(6)  when  £  =  .01,  a  =  .1, 

1,  10,  100  and  Y  =  3.3  and  1,  respectively,  corresponding  to  crack 
propagation  speeds  near  the  glassy  and  equilibrium  values.  G^  is  seen  to 
be  a  montonic  increasing  function  of  the  crack  separation  distance. 

Material  effects,  corresponding  to  the  various  values  of  a,  are  seen  to 
have  a  greater  effect  upon  the  ERR  at  a  large  crack  propagation  speed,  as 
in  Figure  5,  than  at  the  slower  crack  speed  displayed  in  Figure  6.  The 
previous  comments  concerning  6 /a  are  consistent  with  the  behavior  of  G^ 
illustrated  in  these  figures.  Figures  7  and  8  display  G ^  versus  log(6) 
for  varying  crack  speeds  and  a  =  .1  and  100,  respectively.  These 
figures  show  that  the  ERR  is  more  sensitive  to  changes  in  crack  propagation 
speeds  and  crack  separation  distances  at  smaller  values  of  a.  The  results 
of  Figure  7  show  that  except  for  a  small  range  of  slow  crack  speeds,  the 
influence  of  the  two  cracks  upon  the  ERR  decreases  with  increasing  crack 
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propagation  speed.  When  a  is  large,  as  in  Figure  8,  the  reduction  in  the 
0  work  input  into  the  crack  tip  due  to  the  prescence  of  two  cracks  is  only 

overcome  at  large  crack  separation  distances.  Figure  9  illustrates  the 
influence  of  viscoelasticity  in  conjunction  with  crack  separation 
#  distance.  As  might  be  expected,  as  viscous  effects  increase,  G  tends  to 

the  value  of  the  single  crack  problem. 

In  evaluating  (22)  it  was  found  that  the  qualitative  behavior  of 


l/X^Ci/a^)  is  similar  to  that  of  with  the  SIF,  K,  being  similar  to 

the  single  crack  SIF  presented  in  [2] .  For  this  reason  the  normalized  SIF, 


k  =  y-—  ,  which  with  the  results  of  [2]  is  seen  to  be  k  =  — - - - - 

La5  / -  * 

e  e  / 1+3  Xjj(i/a  ) 

2  e 

is  plotted  only  for  the  standard  linear  solid. 


[ 
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Figure  1. 


Figure  2. 


Figure  3. 


Figure  4. 


Figure  5. 


Figure  6. 


Figure  7. 


Figure  Captions 

G2  versus  log(y)  for  y  <  1  and  log((/n+T  -  1 ) / ( /n+T  -  y)) 
for  1  <  y  <  /n+1  for  a  standard  linear  solid  with  6  =  1, 

n  =  10,  e  =  .01,  and  a  =  .1( — ),  i.o( - ),  10( - ), 

100( - •). 

g  versus  log  (y )  for  0  <  y  <  1  and  log((/n  +  1-1 )/(/n +l-y )) 

for  1  <  y  <  /n  +1  for  a  standard  linear  solid  with  6=1, 

n  -  10,  e  =  .01,  and  a  =  .  1(  — ) ,  1.0( - ),  10( - ), 

100( - )  . 

versus  log(y)  for  a  power-law  material  with  6=1, 

n  =  .3,  e  =  .01,  and  a  =  .l( — ),  i.o( - ),  10( - •)  , 

100  ( - ). 

g  versus  log  y  for  a  power-law  material  with  6=1, 

n  =  .3,  t  =  .01,  and  a  =  ,1( — ),  1.0( - ),  10( - ), 

100( - ). 

G2  versus  log(6)  for  a  standard  linear  solid  with  n  =  10 , 

e  =  .01,  y  =  3.3,  and  a  =  .1  ( — )  ,  1.0( - ),  10( - )  , 

100( - ). 

G2  versus  log(6)  for  a  standard  linear  solid  with  n  =  10 , 

e  =  .01,  y  =  1.0,  and  a  =  .1(— ),  1.0( - ),10( - ), 

100( - ). 

G2  versus  log(6)  for  a  standard  linear  solid  with  n  =  10 , 

e  =  .01,  a  =  .1,  and  y  =  .1( — ),  1.0( - ),  10( - ), 

100( - ). 
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Figure  8. 


Figure  9. 


Figure  10. 


versus  log(6)  for  a  standard  linear  solid  with  n  =  10 , 

e  =  .01,  a  =  100,  and  Y  =  .  1  ( — ) ,  1.0( - ) ,  10(—  •— ), 

100( - ). 

G2  versus  log(6)  for  a  standard  linear  solid  with  t  =  .01 

a  =  1,  and  Y  =  1,  and  n  =  .1( — ),  1.0( - ),  10( - ), 

100( - ). 

1/2 

The  normalized  SIF  k  versus  Y/(n+l)  for  a  standard 

linear  solid  with  6=1,  n  »  10,  and  a  -  ,1( — )  ,  1 . 0 ( - 

10. 0( - ). 
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